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1 Introduction 

7— I ■ 

The generation of nonclassical light and its interaction with matter are receiving intense 
Q\ . attention in quantum optics, driven by research opportunities emerging from present control- 

technology regarding atoms and electromagnetic fields. The concomitant, all important details 
of the matter-field interaction have been thoroughly scrutinized, with the Jaynes-Cummings 
model (JCM) playing a pivotal role [U [21 El H]. Despite of the JCM-simplicity, it permits a 
variety of generalizations, applicable to distinct environments and regimes YZ\ . In particular, we 
K ; my nation the work of Buck and Sukunw Q that induced the intensity dependent JCM. 

Because of the commensurability of the Rabi frequencies arising from the model's couplings, 
periodic revivals emerge, absent in the original JCM, with a time-dependent state-vector that 
is periodic itself. As a consequence, any expectation value will share such feature, which leads 
to an enhancement of certain effects that would otherwise be ignored by JCM-practitioners 
[HE]. We wish here revisit the Buck-Sukumar model with information-theory tools so as to be 
in a position to display hopefully interesting details of the model's dynamics. 
These tools have also been the subject of much interest, in particular when they are applied 
in non-thermal setting. In this regard, von Neumann's (NE) [6], linear (LE) [7], and Shan- 
non's entropy (SE) [8j have been frequently used for a variety of quantum systems. It is worth 
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mentioning that the SE involves only the diagonal elements of the density matrix and in some 
cases yields information similar to that obtained from either the NE or LE measures. Other im- 
portant entropic-setting involves semiclassical physics and one employs there the semiclassical, 
phase-space field Wehrl entropy (FWE) [9]. In turn, the FWE has been successfully applied in 
descriptions of different properties of quantum optical fields, such as phase-space uncertainty 
[10] . decoherence (TTJ [12], etc., a theme that will be the focus of our concern in this work. 
Thus, we are led to the concept of Wehrl phase distribution (WPD), that has been extensively 
developed and shown to be a successful indicator of both noise (phase-space uncertainty) and 
phase randomization [13]. Furthermore, the FWE has been fruitfully applied to dynamical 
systems too. In this respect it must mention that the FWE-time evolution in the case of the 
Jaynes-Cummings model has been thoroughly investigated in [T2| (HJ [15] . The FWE i) turns 
out to be more apt for distinguishing amongst states than the NE [13] and ii) is known to yield 
helpful information on atomic inversion processes. Indeed, FWE-studies of a single-trapped ion 
interacting with a laser field with different configurations of the laser field have been consid- 
ered in [16J. We also know now that both (a) the fluctuations of the laser phase and (b) the 
initial-state setting play important roles concerning the evolution of quantifiers like the Husimi 
Q— function, Wehrl's entropy and Wehrl's phase distribution[16j. 

A rather different functional of the probability distribution function (PDF), called Fisher's 
information (FI) [T7] will also be invoked here. FI was originally introduced by Fisher [17] as 
a measure of "intrinsic accuracy" in statistical estimation theory. We will concern ourselves 
in this communication with the Fl-version constructed with the semi-classical Husimi PDF 
[T8| [19j [20]. It has been shown in [21] that FI can be used for evaluating the accuracy limits of 
a quantum measurement because it provides one with meaningful error estimates, even in the 
case of highly nonclassical regimes. This is due to the fact that variances are used to quantify 
the error in quantum measurements (variances and FI are intimately linked via the Cramer-Rao 
bound p2]). The relation between the so-called atomic Fisher information (AFI) and different 
entanglement measures such as von Neumann's, linear, and atomic Wehrl's entropy has been 
investigated in [22J, whose authors found that the entanglement of a two- level atom can be 
measured by using it. Also, FI is used to measure the correlation between the quantized field 
and a Kerr medium [23]. A still new application for FI is found in [21] : it can be employed as 
an information quantifier for the description of the weak field versus strong field dynamics in 
the case of a trapped ion in a laser field. Ref. [23] compared FI, as an information quantifier, 
with von Neumann's and Wehrl's entropies, and provided some analytical Fl-results. 

In the present contribution, that utilizes all the above quantifiers again, our main interest lies in 
a different direction. We wish to investigate the Fl-temporal evolution for a single-qubit system 
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in the presence of an intensity dependent field. Now, given the pertinent PDF we can easily build 
up its marginals and, in turn, a Fisher measure for each of the resulting PDFs. We ask then: can 
all these FIs be used as a quantifiers of the classical correlations and dynamical properties of the 
system ar hand? Additionally, we will focus attention on the effects of i) the laser phase and ii) 
the initial state setting on the evolution of both Wehrl's entropy and Fisher's information. Why 
does this matter? Because these two measures describe (a) interesting semiclassical physics and 
(b) both classical correlations and also quantum entanglement [25j |26| 127] . 
The paper is organized as follows: Section 2 deals with preliminary matters: the basic model 
of a single-qubit in the presence of an intensity dependent field together with the Wehrl entropy 
fundamentals, one the one hand, and the different FIs used here on the other one. Section 4 is 
devoted to the discussion of our numerical results and some conclusions are drawn in Section 
5. 

2 Preliminary materials 
2.1 The model 

We give below the main details for the treatment of a two-level atom interacting with a single- 
mode of the cavity field [5]. The Hamiltonian, in the rotating wave approximation, can be 
written as j5] 

H = UF tfa + ^(|e)(e| - \g)(g\) + Xf (tf\l)(0\ + o|0>(l|), (1) 

where ujf is the field frequency, ua the transition frequency between the upper |0) and lower 
state |1) of the atom, and A the effective coupling constant. The field creation (annihilation) 
operator is a) (a) while fia^a) represents the intensity dependent function of the cavity field 
mode. Restricting ourselves to the functional form f{n) = V a) a, the interaction Hamiltonian 
reads 

^ = A(^|i)(o| + |0)(l|^), (2) 

where ip = aV a) a and $ = V a) act) . The time evolution operator for the effective Hamiltonian 
(2) becomes 



U\t) =exp[-iH T t} = 



3 
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i ; — — ft cos ( T\l ip^ip 



In this last relation T = Xt is the scaled time. The time units are given by the inverse of the 
coupling constant A. We assume (I) that the initial state of the system is the product p AF (0) = 
p A (0) ® p F (0), with our qubit assigned initially to the upper state, i.e., p A (0) = |0) (0|, while 

oo 

(II) the field's initial state is a coherent-one p F (0) = \a)(a\ — ^2 C n (a) C* m (a) \n) (m\ , with 



n,m=0 



C n (a) = a r 



-M 2 



rv. 



The Husimi Q— function Q F of the field-mode, in terms of the diagonal elements of the density 
operator in the coherent-state basis, is 

Qf (M = -Tr^ (P \U(t)p AF (0)tf(t)\p) , (4) 

7T 

where Tr^ means that we trace over the atomic variables. 

Next, we turn our attention to the semiclassical Wehrl-entropy [9] that describes the time 
evolution of a quantum system in phase-space. This entropy, introduced as the classical entropy 
of a quantum state, yields meaningful insights into the dynamics of the system [9] and is defined 
as the coherent-state representation of the density matrix [9l [12] via 

Sw (*) = -/ Qf (P> *) ln Qf *) d2 ^ ( 5 ) 

where d 2 j3 = \(3\d\f3 \ dQ . We can specialize things by recourse to the Wehrl phase distribution 
(Wehrl PD), defined to be the phase density of the Wehrl entropy [13l[T6], i.e., 



Se(t) = - J Q F (f3,t)\nQ F (f3,t)\f3\d\f3\ (6) 

where 6 = arg (/?). 

2.2 Fisher Information 

The Fisher information measure (FIM) for any PDF /(x) can be cast in the fashion [281 129] 
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and is encountered in many physical applications (see, for instance, [30]- jH], and references 
therein). The FIM associated to Husimi distributions Qp(Xi, X 2 ,t) is defined as [4"5] 

I poo roo 

I F (t) = J Q F (X 1 ,X 2 ,t)T(X 1 ,X 2 ,t)dX 1 dX 2 , 

Q F ((3,t)T((3,t)d 2 f3, (8) 



where T(X 1 ,X 2 , t) and T (/3,t) can be written in terms of the phase space parameters, yielding 



r(x 1 ^,^t(^(*) ^g 1 '^ t)) )', 



j=l k=l 
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cos I H — 



ain(g (M) 



[k — 2 + /3 (k — 1)] d\p (k — 2) + Q (k — 1)] 



V 



/ 



with 



and 



*x i (t) = y/(X j (t)*)-((X J (t))) 2 , 



(9) 



(10) 



/oo /*oo 
/ X j Q(X 1 ,X 2 ,t)dX 1 dX 2 . 
■oo J —oo 

We also consider, as a dynamical measure, the quantity 

POO 

Ie(t) = - / Q H (Mr(/M)|0|d|0|. 



(12) 



It is worth noting that the definition (jSj) is given in analogy to that of the field Wehrl entropy 
(the special case t = being lp(0) = 2) so that the corresponding Fisher's phase distribution 
can be cast, in terms of the error function erf (x) = e~ y2 dy, as 

1 



Io(0) = — exp( 



2 2\ 

j -a j 



x 



{xa/tt[1 + erf (x)]/i + exp (-x 2 ) / 2 }, 



where x = a cos(9) and 



fj = a 2 -x 2 + J -, j= 1.2. 



(13) 



(14) 



In correlations terms the bipartite system becomes uncorrelated whenever Jf(0) ~ 2, this value 
representing the lower bound for Ip. One has 



I F (0) = S w (0) + 1-1H7T, (15) 

and 

J e (0) = S e (0) - ^ ex P ("« 2 ) |l + xv / ^(l+erf(x))exp(x 2 )|. (16) 

Equations (|15p - ffTBI) establish the connection between Fisher's information measure (FIM) 
and Wehrl's entropy at T — 0. Notice that at this time Ip — Iq = 1 — In tt = constant, which is 
a counterintuitive result, since one expects their sum to be approximately constant |44] . This 
curious result is due to the periodicity of the evolution. 
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3 Results 



We start now the presentation of our numerical results. We will see that the coherent state 
parameter a, representing the square root of the mean-photon number, greatly influences the 
dynamics, as can be clearly appreciated in Figs. [H that depict, respectively, If and Sw as a 
function of T and a [(a) and (b)] together with their projections on the a — T— plane [(c) and 
(d)] (T is a "scaled" time). Both the inherent periodicity of the dynamics and the long living 
correlation between the single qubit and the coherent field are clearly visible. They increase as 
the phonon-number grows. Both quantifiers exhibit the periodicity of the system. 
In figure [2] we plot the FI- and Wehrl- time evolutions together with the associated Y — variance. 
For typographical simplicity, we set in the graphs Y = X 2 . One chooses three values of the 
a— parameter, namely, = 1,2,3, respectively. In order to ensure good accuracy, the behav- 
ior of the Fisher information Ip(t) has been determined using an appropriate scale so as to 
meaningfully compare it to Wehrl's entropy. FI's behavior is clearly dominated by the variance 
component a 2 X2 {t) [Cf. Eq. [TT] . 

Fig. [3] illustrates the If vs. Sw behavior, which is surprising indeed. The second is a global 
measure, while the former is a local one [IB]. One expects them to behave in opposite fashion 
[48J. Here though, a monotonous comportment is apparent, for the first time as far as we know. 
Rapid fluctuations, combined with strong derealization are thus generating a bizarre scenario. 
There is more to be said, though. Wherl's entropy was invented as a measure of derealization 
in phase-space. For the minimum possible mean photon number (corresponding to a = 1), 
quantum indeterminacy is maximal and both quantifiers grow together. We may understand 
this behavior if we realize that, as indeterminacy augments (and so does Sw) fluctuations rise, 
forcing Ip to increase as well. Enters now a bit of rather interesting physics. The route to 
classicality is paved by i) a growing mean photon-number and ii) a stabilization of Wehrl's 
entropy since derealization stops augmenting. This scenario begins to insinuate itself at a = 2 
and becomes fully installed already at a = 3! Thus, Sw cannot grow, but nothing impedes If 
to increase. What we are really watching in Fig. [3]-(c) is the emergence of the classical limit. 
The above considerations receive a boost via Fig. 0], that depicts the FIM- Wehrl behavior 
versus the mean phonon number. In a sense, things return to "normality", as the quantifiers 
behave differently here, as expected. Ip always grows with a entailing that, as one intuitively 
understands, errors diminish as particle-number grows. Wherl's measure has a peak at about 
a = 3. This can also be understood. Sw measures our ignorance about localization in phase- 
space, which, as it should, diminishes as a grows. 
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3.1 Cramer-Rao bound 



The "true" Fl-informational content is conveyed by the Cramer- Rao inequality (CR). Indeed, 
this is its most important property. We recapitulate in one-dimension, for simplicity's sake. If 
the classical Fisher information associated with translations of a one- dimensional observable x 
with corresponding probability density f(x) is [2E1 HS| 

din N 2 

then it obeys the above referred inequality, namely 



/.= / dxf(x) (— ^) , (17) 



(Ax) 2 > I' 1 

involving the variance of the stochastic variable x 



> 2 

„2 



(Ax) 2 = (x 2 ) - (x) 2 = J dxf(x)x 2 - dxp(x)xj . (19) 

We remark that the derivative operator significantly influences the contribution of minute local 
/—variations to FI's value, so that the quantifier is called a "local" one. Note that Wehrl's 
entropy decreases with skewed distributions, while Fisher's information increases in such a 
case. Local sensitivity is useful in scenarios whose description necessitates appeal to a notion of 
"order" . For our present purposes we deal with a time-dependent CR, that, in self-explanatory 
notation reads 



where 

A? = (/? 2 > - (/3) 2 

(/3 s ) = C io°° Qf t) \P\ d \P\ dQ -> s = 1 ^ 

and we display in Fig. [5] some results. We see that the CR product oscillates with time and 
rapidly decreases as a grows. At a = 3 the product almost reaches its lower bound of unity, 
reconfirming the fact that at that a— value classicality has been reached. 



4 Conclusions 

We have here considered, from an information-theory viewpoint, the dynamics of a single- 
qubit system. Our information-quantifiers were the Wehrl entropy, a phase-space derealization 
measure, and Fisher's information. These two quantities aptly illustrate on the complicated 
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dynamics at hand. The main characteristics of the problem are governed by the mean phonon 
number and the intensity dependent field. 

We have performed an extensive numerical analysis, illustrated via a variety of graphs. Peri- 
odicity is a main feature, specially as the mean-phonon number grows. Long living correlations 
between the qubit system and and coherent field are clearly appreciated. The monotonous 
growing of the Fisher measure as the Wehrl entropy grows is a counterintuitive feature that 
we have detected. This is a surprising facet because it is well known that whenever Fisher 
grows, Shannon or Wehrl entropies decrease J28JIH]. Notice that FI measures gradient content 



\28\ 144] while Wehrl's measure is a derealization- indicator [44J. However, the physics of the 
phenomenon can be understood, as explained in the preceding Section. Rather unexpectedly, 
we get illuminating insights into the emergence of classicality, in line with very recent findings 
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Figure 2: Time volution of (a) Fisher' information Jp(t), (b) Wehrl's entropy Sw(t), and (c) 
variances a\ At). We set Y = X 2 . We also plot curves for different values of the square root of 
the mean photon number a, namely, dotted curve for a = 1, dashed curve for a = 2, and solid 
curve for a = 3. 
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Figure 3: Illustration of the Ip vs. Sw behavior for parameters values for different values of 
the square root of a, where figure (a) a = 1, figure (b) a = 2, and a = 3 for figure (c). 
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Figure 4: We depict the FIM-Wehrl behaviors versus the square root of the mean photon 
number. 
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Figure 5: The Cramer- Rao product Ip{t)^ as a function of the scaled time T for different 
values of a where the dotted curve a — 1, dashed curve a = 2 and the solid curve a = 3. 
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